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Abstract 



It is extraordinary that a number of observations indicate that we hve in 
a spatially flat, low matter density Universe, which is currently undergoing a 
period of accelerating expansion. The effort to explain this current state has 
focused attention on cosmological models in which the dominant component 
of the cosmic energy density has negative pressure, with an equation of state 
w > —1. Remarking that most observations are consistent with models right 
up to the w = —1 ox cosmological constant (A) limit, it is natural to ask 
what lies on the other side, at w < —1. In this regard, we construct a 
toy model of a "phantom" energy component which possesses an equation 
of state w < —1. Such a component is found to be compatible with most 
classical tests of cosmology based on current data, including the recent type 
la SNe data as well as the cosmic microwave background anisotropy and 
mass power spectrum. If the future observations continue to allow < — 1, 
then barring unanticipated systematic effects, the dominant component of the 
cosmic energy density may be stranger than anything expected. 
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Arguments have been put forward that we hve in a spatially flat, low matter density 
Universe which is currently undergoing a period of accelerating expansion. If the observa- 
tional evidence upon which these claims are based are reinforced and strengthened by future 
experiments, the implications for cosmology will be incredible. It would then appear that 
the cosmological fluid is dominated by some sort of fantastic energy density, which has neg- 
ative pressure, and has just begun to play an important role today. No convincing theory 
has yet been constructed to explain this state of affairs, although cosmological models based 
on a dark energy component, such as the cosmological constant (A) or quintessence (Q), 
are leading candidates. At this stage we are lead to notice that the parameterization of the 
dominant energy component as a fluid with an equation of state (defined as the ratio of 
pressure to energy density) w{= p/p) > —1 leads to the curious situation that most obser- 
vational constraints are consistent with models that go right up to the w = —1 border [|l|. 
It is natural to ask what lies on the other side of this boundary. The focus of this paper is 
the investigation of cosmological models with a dominant component for which w < —1. 

We begin by constructing a classical cosmology: a spatially flat FRW space-time filled by 
a cold dark matter (CDM) component and a "phantom" (P) energy component. (A phantom 
is something which is apparent to the sight or other senses but has no corporeal existence - 
an appropriate description for a form of energy necessarily described by unorthodox physics, 
as will be seen later.) The phantom energy has positive energy density, pp > 0, but negative 
pressure, such that Pp+Pp < 0. (It is no fatal flaw for a component to violate the dominant 
energy condition for a finite time, as can arise from a bulk viscous stress due to particle 
production 0.) It immediately follows from the equation of state w < —1 that the phantom 
energy density grows with time. If this energy has begun to dominate today, it must have 
come on like a bolt from the blue. Taking the constant value w = —3 as an example, then 
Pp oc where a is the expansion scale factor. So if Qm ~ 0.3 today, then the Universe 
contained 90% CDM at z ~ 0.4 as opposed to 2 ~ 1.8 in a A-dominated cosmology. Hence, 
the phantom energy exerts its influence very late. We now turn to quantify the remarkable 
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effect on the cosmology of a component with super- negative equation of state, w < —1. 

The expansion scale factor grows rapidly when the phantom energy comes to dominate 
the Universe, as a/a actually grows and the deceleration parameter qo = (1 + 3wQp)/2 
becomes very negative. In fact, the scale factor diverges in finite time: if the expansion 
is matter dominated until the time tm, then we can write the phantom energy-driven scale 
factor as a(t) = a{tm)[—w + {l+w)t/tm]'^^'^''^^'^^ at t > tm- For w = —3 then a diverges when 
if: = 3tm/2. Unless Qm ^ 1, however, to < tmw/^l + w) and so the limiting cosmological 
time occurs well after the present day. 

The expansion age, and similarly the horizon distance in a phantom energy cosmology 
are larger than in the analogous A model. In Figure we show the age for a sequence of 
w < —1 models. In the range of interest, near Q^n ~ 0.3 — 0.4, we see that the phantom 
energy can increase the age by up to ~ 30% over the w = —1 age. For a given value of Qrn, 
as w becomes very negative the phantom energy becomes important later and later with 
diminishing effect on the age. In fact, for = 0.3, the bound on the age in units of the 
Hubble time is -f^o'^o < 1-2. 

The volume - red shift relationship is demonstrated in Figure 0. For the same matter 
density, the phantom energy model gives the largest differential number of objects per red 
shift interval. Although evolutionary effects are important in cosmological tests based on 
this relationship, if all other features are held fixed, the phantom model will predict more 
strong gravitationally lensed quasars than Q or A models. 

The magnitude - red shift relationship is demonstrated in Figure ^. The predictions for 
several cosmological models have been shown along with a summary of the observational 
results. All other parameters being equal, one expects high red shift supernovae to be 
dimmer in a phantom energy Universe. 

Considering the constraint to the allowed range of w and flm due to the supernovae, 
we find that phantom energy models even with very negative w are in accord with the 
observations. We see in Figure ^ that the contour region is extensive for w < —1, prefering 
a slightly higher matter density than in Quintessence models. An alternative way to visualize 
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this constraint on phantom models is to construct a parameter space in terms of the red 
shift at which the energy density drops to 90% matter, signalhng the end of the matter- 
dominated epoch, rather than the constant equation of state w. The right-hand panel above 
makes clear that an important difference with other dark energy models is the very late 
onset of the phantom energy. 

The next challenge is to construct a microphysical model of the phantom energy in 
order to consistently determine the impact on the cold dark matter, baryon, and photon 
fluctuation spectra. 

We describe the system of gravitation, normal matter, and phantom energy with the 
Lagrangian L = —R/lGrcG + Lfn + Lp where represents the Lagrangian for the "normal" 
matter, CDM plus the standard model particles and fields. Since it is not possible to achieve 
w < —1 with a free scalar field, and we wish to avoid interactions with other fields in order to 
keep the energy "dark," we start with the unorthodox Lagrangian Lp = — 9^09'^0/2 — V{(j)) 

(with metric signature H ). The important point is that we have switched the sign 

of the kinetic term in the scalar field Lagrangian Lp. In this way, Pp = —0^/2 + V and 
Pp = —0^/2 — so that w < — 1 is attained. The equations of motion are + 3H(j) = +V(j,, 
so that the field will tend to run up, not down, a potential towards larger energy. The 
key ingredient for this crude toy model, and the more detailed model we introduce later, is 
the non-canonical kinetic energy term, examples of which occur in super gravities and in 
higher derivative theories of gravity 

The cosmological spectrum of fluctuations in the phantom field develop in a fashion 
similar to the case of Quintessence 0. The Fourier transform of the linearized perturbation 
equation is 6(j) + 3H6(f) + {k"^ — V^,^^)^^ = —h(j)/2, where h is the synchronous gauge metric 
perturbation. In this equation, the sign of T^^^^ is different than in the standard case. Just 
as for Quintessence, perturbations on small scales, for k"^ ^ |^(/)</>|) are suppressed. On 
larger scales, however, the phantom energy develops inhomogeneities in response to the 
surrounding matter and density perturbations. If the effective mass, (/c^ — V^^Y^"^, should 
become imaginary then 50 would develop a growing (tachyon) instability. However, if we 
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implement a constant equation of state, then V^ip = |(1 — w)[H — + w)], which is 

negative for a hmited range of values of w and flm- Consequently, the effective mass is real, 
and there is no instability for constant w. 

A more versatile phantom model can be constructed using a Lagrangian with a non- 
standard dependence on the field gradients, Lp = a(0)(V0)^ + 6(0)(V0)'^ + motivated by 
the appearance of higher derivatives in string theory. This approach has been employed in 
the k-inflation and k-essence scenarios and it has been demonstrated that such 

models can also be used to describe an energy component with a super-negative equation of 
state, Wp < —1, which is also stable to perturbations (no tachyonic instability). Details of 
the treatment of cosmological perturbations in the k-essence scenario are given in the 
approach is sufficiently general that the case of a phantom is simply accommodated. Hence, 
starting from a theory giving the functional dependence of Lp on and (V0)^, one can 
readily determine the stress-energy p{z), p{z) and the sound speed of perturbations Cp{z) as 
they evolve in time or red shift. 

Turning this around, one may instead dictate the histories Wp{z) and Cp{z) in order to 
specify the phantom dark energy model. In the present investigation we focus on constant- 
Wp and Cp = 1 models. Although this is not a generic prediction of the dynamical models 
described above, the observational differences are not expected to be very large (as some 
exploratory analyses have confirmed). Since the observational effects rely on a substantial 
fractional phantom energy density, and since p{z) is actually growing with time, then the 
relevant values of Wp{z) and Cp{z) are the mean values between the redshift when the phantom 
began to dominate and the present. For Wp < —1, this is a narrow range in redshift, leaving 
little opportunity for the effects of time- variation to be distinct. As we have modified 



CMBfast |TT| to include a phantom dark energy component based on both approaches, we 
now proceed to discuss the resulting perturbation spectra. 

The main feature that distinguishes the phantom energy case from A or Q is that the 
onset of phantom energy dominance happens at the very last moment - so late that most 
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evolutionary effects which occur in A and Q models are suppressed. In the cases illustrated 
in Figure |^ the strength of the late-time ISW effect diminishes as w becomes more and more 
negative, because the expansion is CDM-dominated until later and later. On small angular 
scales the locations of the acoustic peaks shift to higher multipole moments as w becomes 
more negative, due to the increased distance to the last scattering surface. This may be 
an important factor in attempts to develop a phantom dark energy cosmological model to 
accommodate the newer CMB results, also represented in Figure |^. A more detailed analysis 
will be forthcoming. 

The phantom-dominated background evolution has an important effect on the fluctuation 
spectra, as well. Making the assumptions that the matter component of the cosmological 
fluid carries a spectrum of scale-invariant perturbations generated by inflation, and that, 
similar to Quintessence, the phantom energy itself does not fluctuate on scales well below 
the Hubble horizon, then we may follow the growth of linear perturbations in the CDM and 
baryons. The growth suppression factor g = D{a)/a is shown in Figure ^, where we see 
that perturbations grow as D oc a until very late, owing to the very late time at which the 
phantom energy begins to dominate. Hence, the evolution is very similar to the standard 
CDM scenario, but with a lower matter density. 

Next, we consider the behavior of the mass power spectrum for w < —1. Because the 
phantom field does not fluctuate on scales well below the horizon, the overall shape of the 
mass power spectrum is well-described by the BBKS parameterization (See for 



instructions on how to adapt the fitting functions to quintessence). The normalization is 
obtained by comparison with the CMB. Using the amplitude of the first acoustic peak, then 
the normalization is similar to a ACDM model with the same cosmological parameters, but 
with a CTg higher by a ratio of the growth suppression factors Qp/qa. Since QAi^m = 0.35, a = 
1) = 0.8, and gp is at most unity in the limit of w ^ —1, then ag is higher by up to 20%. 
Using COBE for the normalization, on the other hand, the effect of the late ISW and the 
direct fluctuations in the phantom field must also be taken into account. However, owing to 
the late onset of the phantom, there is very little late ISW and the direct fluctuations are 



also negligible. As mentioned earlier, the low-i spectrum is very flat, as in SCDM, and up to 
~ 5% lower than in the comparable ACDM model. Hence, the COBE normalized spectrum 
will result in a power spectrum with a as up to ~ 20 + 5% higher than for A. 

In Figure |^ we show the prediction and observational constraint on ag for a sequence of 
models varying in w. While the amplitude of the CMB-normalized mass power spectrum 
grows slightly as w becomes more negative, the implied as based on the observed abundance 
of clusters grows as well. (The cluster abundance calculation is based on an extension of 



IJ] to the case oi w < —1.) We see that the intersection of the 2a regions of as grows 
with decreasing w. Further constraints based on large scale structure one may consider, 
but which are beyond the scope of this paper, include the evolution of the abundance of 
rich clusters and the growth rate of the linear fluctuation spectrum. Our lesson here is that 
phantom dark energy models do exist which, upon first inspection, satisfy constraints based 
on the CMB and the mass power spectrum. 

In summary, we have investigated the properties of cosmological models in which the 
dominant energy density component today has an equation of state w < —1. We have 
demonstrated the impact on the cosmological age, the volume - red shift and magnitude - 
red shift relations, the CMB, and the mass power spectrum, finding broad agreement with 
current observational constraints. Current data appears to be consistent with a phantom, 
although a more careful analysis will be necessary to quantify the observational status of 
w < -1. 

While there is no single, best method for determining if the equation of state is strongly 
negative, combinations of measurements of the type described in this paper will substantially 
improve our understanding of the dark energy. One of our central points is the importance 
that such analyses are open to the possibility of a phantom dark energy - unjustified biases 
and priors can lead to a gross misinterpretation of the observational evidence. (See for 
an analysis of the pitfalls of assuming a constant w or w > —1.) 

Of course, we do not want to overlook the distinct possibility that the observational 
evidence in favor of such a strongly negative equation of state is simply a phantom - that the 



apparent accelerating expansion is due to more conventional, though unanticipated causes 
{e.g. dust or evolution for the SNe). If these systematic effects can be eliminated, and 
the data continue to support w < —1 then the implications for fundamental physics would 
be astounding, since w < —1 cannot be achieved with Einstein gravity and a canonical 
Lagrangian. It is premature to shift attention towards building w < —1 models, but it is 
important to be aware of the properties and implications of models in each direction of the 
cosmological parameter space. 

As has been discussed elsewhere, current observations suggest the presence of a dark 
energy component with w ^ —0.5 IQ, with many constraints pushing towards w = —1: the 
observations are teetering at the edge of a previously unfamiliar boundary. 
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red shift: z 

FIG. 2. The volume - red shift relationship is shown for phantom energy models with 
w = -3, -3/2, ACDM with w = -1, QCDM with w = -1/2, and CDM. All the dark energy 
models have il^ = 0.3. 
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FIG. 3. The magnitude - red shift relationship is shown for the redshift-binned type la SNe data 
(HZT- 0; SCP- [H; binning- [||]), alongside the predictions of various cosmological models. 
Both the phantoms {w = — 3/2, flrn = 0.4; —3/2, 0.3) and the A model (—1, 0.3) provide good fits 
to the data (low x^/d.o.f.). The magnitudes are calculated relative to an empty, open Universe. 
The light, dashed lines are for pure phantom, deSitter, Milne, and Einstein-deSitter, from top to 
bottom. 
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equation of state: w red shift: z 

FIG. 4. The contraint on the allowed values of 0,m in phantom and quintessence dark energy 
models is shown as a function of w, or alternatively the red shift at which matter-domination ends, 
when = 0.9. We have traversed the w = —1 boundary, finding that there are phantom energy 
models in accord with the observations. 
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multipole / 

FIG. 5. The CMB anisotropy spectra are shown for various cosmological models with equal 
input power, to demonstrate the effect of w < —1. On large angular scales, the strength of the 
late-time ISW effect diminishes as w becomes more and more negative, because the expansion is 
CDM-dominated until later and later. The locations of the acoustic peaks shift to smaller angular 
scales as w becomes more negative, due to the increased distance to the last scattering surface. The 
models shown all have $7^ = 0.35, il^/i^ = 0.02, and h = 0.70. The horizontal axis is logarithmic 
for 2 < / < 200 in order to give enough space to both large and small angular scale features. 
Comparing the curves with a compilation of current CMB data |2^ suggest a more negative dark 
energy equation-of-state may allow for a better fit to the small angular scale data. 
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FIG. 6. The amplitude of the growth suppression factor g = D{a)/a versus the scale factor is 
shown for various cosmological models. We see that the perturbation growth is slowed later, with 
weaker effect, in the phantom energy models. The models shown all have Vtm = 0.3. For SCDM, 
D{a)/a = 1. 
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FIG. 7. The cluster abundance constraint is shown for a sequence ofQ {w > —1), A {w = — 1), 
and phantom (w < —1) models with $7^ = 0.3, h = 0.65, Qfyh'^ = 0.02, n = 1. The upper shaded 
region shows the 2a range of fjg based on the observed abundance of x-ray clusters, as interpreted 
for a range of w (based on an extension of the work in Ref. The lower shaded region shows 

the predicted fig for COBE- normalized models. The intersection of the two regions grows as w 
becomes more negative. The letters P and Q indicate the Phantom and Quintessence regions, to 
the left and right of w = —1. 
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